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Motivated by experiments on ultracold atoms which have realized the Haldane model for a Chern 
insulator, we consider its strongly correlated Mott limit with spin-1/2 fermions. We find that slave 
rotor mean field theory yields gapped or gapless chiral spin liquid Mott insulators. To study com¬ 
peting magnetic orders, we consider the strong coupling effective spin Hamiltonian which includes 
chiral three-spin exchange. We obtain its classical phase diagram, uncovering various chiral magnetic 
orders including tetrahedral, cone, and noncoplanar spiral states which can compete with putative 
chiral quantum spin liquids. We study the effect of thermal fluctuations on these states, identifying 
crossovers in the spin chirality, and phase transitions associated with lattice symmetry breaking. 
We also discuss analogous effective spin Hamiltonians for correlated spin-1/2 bosons. Finally, we 
point out possible experimental implications of our results for cold atom experiments. 


I. INTRODUCTION 

Momentum space topology is key to our understand¬ 
ing of phases such as topological insulator^II^ or quantum 
anomalous Hall insulatorsThe interplay of momen¬ 
tum space topology and local real space interactions is 
expected to lead to a rich variety of correlated insulating 
phasesThis has motivated an extensive investiga¬ 
tion of the time-reversal invariant Kane-Mele-Hubbard 
model or more realistic variantwhich provide the 
simplest examples of interacting quantum spin Hall insu¬ 
lators. Recently, cold atoms in a shaken optic al lattic^^ 
or laser-induced tunneling or transitions^^^^ have been 
proposed for realizing topological states of matter, and 
employed to realize the Haldane honeycomb model of a 
Chern insulator (Cl) and Chern bands in the square 
lattice Hofstadter model with flux 7r/2 per plaquette.^^ 
Preliminary experiments suggest that interactions in the 
Haldane model may not lead to excessive heating.l^ Mo¬ 
tivated by this, we explore strong correlation effects in 
the Haldane model for spin-1/2 fermions and bosons. 

The Haldane model is defined on the two-dimensional 
(2D) honeycomb lattice in Fig. 1(a), with a real near¬ 
est neighbor hopping amplitude ti, and a complex 
next-neighbor hopping t 2 C^ which breaks time-reversal 
symmetry.l^ For ^2 = 0, the energy dispersion is identical 
to graphene, supporting two inequivalent massless Dirac 
fermion modes. For small ^2 7^0, and 0 7^ 0,7r, the Dirac 
fermions acquire a mass, leading to a quantized Hall ef¬ 
fect (Jxy = jh at half-filling. For spin-1/2 fermions, 
with each species at half-filling, jh. What 

is the fate of this quantum Hall insulator when interac¬ 
tions lead to Mott localization? Does the Mott insulator 
support unusual magnetism or spin liquid ground states? 

For spinless fermions, nearest-neighbor repulsion 
induces a to pologically trivial charge density wave 
insulatorj ^^ l ^^ l while spinless bosons with a local Hub¬ 
bard repulsion form a plaquette Mott insulator with 
loop currents.l^ A slave-spin approach to the Haldane- 


Hubbard model suggests the appearance of fractional¬ 
ized Chern insulators^ or topological Neel order at mod¬ 
erate coupling with a staggered sublattice potential.!^ 
An alternative slave rotor mean field theory of spin-1/2 
fermions, in which Hubbard repulsion localizes only the 
charge degree of freedom, led to a me an fie ld descrip¬ 
tion of a gapped chiral spin liquid (CSL).l^Sinil Qutzwiller 
projected wave function studies show that such a CSL 
would b e in th e same phase as the u = l/2 bosonic Laugh- 
lin state.^^^^ More generally, such CSLs have al so be en 
shown to be induced by chiral 3-spin interactionsor 
appear naturally for SU(N) fermion^^ at large N. 

However, previous work has not considered the var¬ 
ious types of magnetic orders which can compete with 
quantum spin liquid phases. Indeed, even the type of 
short-range spin correlations in the correlated Mott in¬ 
sulator, which could potentially be probed in cold atom 
experiments, have not been explored. Finally, since the 
topological order associated with such chiral spin liquids 
will not persist at nonzero temperature, it is important to 
gain an understanding of the thermal fluctuation effects 
on the Mott insulating state if we are to make connec¬ 
tion with experiments. However, previous studies of this 
model have all been restricted to the T = 0 limit. 

This motivates us to revisit the effect of strong corre¬ 
lations on the Haldane-Hubbard model. Our key results 
are the following, (i) We extend the slave rotor mean field 
theory of the fermionic Haldane-Hubbard model to larger 
^2- We show that this parton construction leads to CSLs 
with gapped or gapless bulk spinon excitations which de¬ 
scend from Cl or Chern metal (CM) phases. The gap¬ 
less spin liquids display multiple spinon Fermi pockets. 
Both types of spin liquids are expected to support gap¬ 
less spinon edge states. Realizing the Haldane-Hubbard 
model with a broad range of ^2 would thus allow one 
to explore the physics of gapped as well as gapless chi¬ 
ral spin liquid phases, (ii) The gapped or gapless spin 
liquids uncovered in the slave rotor parton description 
may be unstable to magnetic ordering due to fluctua- 
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tion effects beyond mean field theory. To study mag¬ 
netic instabilities in the strongly correlated limit, we con¬ 
sider the strong coupling limit of the Haldane-Hubbard 
model, and derive and study the effective spin Hamilto¬ 
nian which has chiral 3 -spin interactions. This leads us 
to uncover a rich variety of competing chiral magnetic 
states such as tetrahedral, cone, and noncoplanar spirals 
which could compete with spin liquid phases. Quantum 
melting such phases will typically lead to broken sym¬ 
metry states. However, unlike other states, we find that 
the tetrahedral state has spin correlations which respect 
all lattice symmetries and a large uniform chirality. We 
tentatively identify this state as the classical parent state 
of a chiral quantum spin liquid. Such tetrahedral states 
have also been recently uncovered in a honeycomb lat¬ 
tice Hubbard model at intermediate correlation^^, (hi) 
Turning to physics at nonzero temperature, we discuss 
thermal fluctuation effects and lattice symmetry break¬ 
ing thermal phase transitions associated with many of 
these states using a combination of analytic Landau the¬ 
ory arguments supplemented by classical Monte Carlo 
simulations. Such thermal transitions may be observable 
in experiments, (iv) Finally, we highlight some analogous 
results for spin-1/2 bosons which can also be studied in 
cold atom experiments, but which have not been theo¬ 
retically explored. 


II. FERMI HALDANE-HUBBARD MODEL 

For fermions with spin-1/2, the Haldane Hubbard 
model is described by the Hamiltonian 

Hnn = iclc^a + h.c.)-t2j2 (e^^'^^clcj^+h.c.) 

rii^riii ( 1 ) 

i 

where (.) and ((.)) denote, respectively, first and second 
neighbors, and Uij = ±1, depending on whether we hop 
along or opposite to the arrows shown in Fig. 1 (a), and 
U is the local Hubbard repulsion. 

For 1 / = 0 , diagonalizing the Hamiltonian leads to 
Chern bands with Chern numbers (7 = ± 1 . Depending on 
12/ti and 0, these bands may overlap in energy, leading 
to a CM with a non-quantized Hall effect, or be well sep¬ 
arated in energy leading to a Cl (equivalently, a quantum 
Hall insulator) with a Hall conductivity jh per 

spin. This phase diagram is shown in Fig. 1 (b), along 
with (Jxy which is quantized in the Cl but non-quantized 
in the CM. (Note that the ^2 = 0 axis is a Dirac semimetal; 
on this singular line, the fiux plays no role.) 

HI. MEAN FIELD CHIRAL SPIN LIQUIDS 

Upon increasing the Hubbard repulsion, this topolog¬ 
ical band metal or band insulator will transition into a 


t/ = o 




FIG. 1: (a) Honeycomb lattice showing the two sublattices 
with hoppings ti,t 2 and phase 0 defining the noninteracting 
Haldane model, with large and small triangles enclosing fluxes 
±30 and —0 respectively, (b) Phase diagram at half-filling for 
each spin species showing Chern band metal and Chern band 
insulator, and the Hall conductance per spin axy in units of 
jh. 


correlated Mott insulator once interactions exceed the 
bandwidth. The critical repulsion Uc = 877|A|, where 
K = A(ti,t2,0) is the kinetic energy per site per spin 
of the noninteracting bands at half-filling. 

Slave rotor mean field theor}^^, provides a concrete 
realization of the Mott phase and the Mott transition; 
it is a parton construction in which the spin and charge 
degrees of freedom are carried by two partons, a neutral 
spin- 1/2 fermion and a charged spinless rotor. The local¬ 
ization of the rotor angular momentum corresponds to 
the transition into the Mott ins ulato r phase. A single¬ 
site slave rotor mean field theor3J^^^ yields 77 = 1 . How¬ 
ever, a comparison with quantum Mont e Car lo studies 
of the honeycomb lattice Hubbard model^^^ at ^2 = 0 , 
and variational numerical studies of the triangular lat¬ 
tice Hubbard modeP^ with ti = 0 = 0 show that the 
single site slave rotor theory overestimates 77; the Mott 
transition occurs at a smaller renormalized 77^0.6-0. 7 . 

Fig. shows the phase boundaries which delineate the 
topological CM or Cl from the correlated Mott insula¬ 
tor (assuming Uc = 877IA| with 77 = 0 . 65 ) at an interac¬ 
tion strength U/ti = 10 . Within slave rotor mean field 
theory, the spinon band structure in the Mott insulator 
resembles the Cl or CM band structure from which the 
Mott insulator descends. Thus, depending on ^2,0, one 
could have gapped or gapless spinon excitations in the 
bulk, accompanied by gapless chiral edge modes; we in¬ 
dicate these mean field states as ‘gapped CSL’ and ‘gap¬ 
less CSL’. The spinon band structure mean field theory is 
simply inherited from the parent noninteracting fermion 
dispersion. Since the parent Chern metal phase exhibits 
multiple Fermi pockets, as shown in Fig. the gapless 
CSL will inherit these as spinon Fermi pockets. Gapless 
chiral spin liquids have also been recently discussed in 
kagome antiferromagnets.E^ 

Going beyond mean field theory, we expect the ‘gapped 
CSL’ to become a genuine spin liquid with semion excita- 
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FIG. 2: Slave rotor mean field theory of the Haldane-Hubbard 
model showing the topological Chern metal (CM, green), 
Chern insulator (Cl, brown), and topological spin liquid Mott 
insulators (black) at Vjh — 10. The topological Mott insula¬ 
tor could have gapless or gapped spinon excitations depend¬ 
ing on whether it descends from a CM or Cl phase, leading 
to gapless or gapped chiral spin liquids. An example of the 
parent Fermi surface of the non-interacting model, which is 
inherited by the spinons at mean-field level, is shown at a 
particular point. The slave rotor result has been rescaled by 
r] = 0.65 (see text). 


tions and topological degeneracy) ^^ ! ^^ * While these CSLs 
are expected to be stable for weak gauge fluctuations, 
strong gauge fluctuations might drive magnetic ordering 
instabilities, or, spi non p airing instabilities in the case 
of the gapless CSL.E3SI Since slave rotor theory does 
not capture competing magnetic ordering tendencies of 
the Mott insulator, we address such potential competing 
phases using a strong coupling expansion.!^ 


IV. COMPETING MAGNETIC ORDERS 

At large I/, a standard derivatiorP^^^ leads to an ef¬ 
fective spin Hamiltonian 






U 

i^j) 

2Atlt2 


u 


{{ij}) 


24:th 


XAsin$A--^E XAsin$A( 2 ) 


small—A 


big-A 


where xa = • Sj x S/. is the scalar spin chirality oper¬ 

ator defined around triangular plaquettes, with the sites 
{ijk} being labelled going anticlockwise around the tri¬ 
angles. being the flux enclosed by the correspond¬ 
ing triangle. As shown in Fig. 1(a), the set of triangles 
includes small triangles that enclose flux ^a = —0 , big 
triangles within a hexagon which enclose a flux $a = 30, 
and big triangles defined around a site which enclose a 
flux ^A = —30. 

Treating the spin as a classical vector, we have ob¬ 
tained the magnetically ordered ground states of i^spin 
using numerical simulated annealing and variational spin 


configurations. This approximation may be formally set 
up by scaling the 2-spin exchange couplings by 1/(45'^) 
and the 3-spin exchange couplings by 1/(85'^) and then 
taking S ^ oo which leads to a classical spin model with 
quantum fluctuations suppressed by 0{1/S). Below we 
will use this classical approximation for S' = 1/2 as an 
uncontrolled approximation, which is nevertheless known 
to work well in many cases in predicting the correct mag¬ 
netic order in 2D quantum spin models. For plotting the 
ground state spin configurations and structure factors, we 
view the honeycomb lattice as a square lattice with 1/4 
deleted bonds (‘brickwalF lattice), similar to the geome¬ 
try realized in optical lattice experiments.^^The resulting 
phase diagram, shown in Fig. 3 for Ujt = 10, contains 
six phases: Neel, Tetrahedral, Cantellated Tetrahedral, 
Spiral, Cone-I and Cone-II, with a line indicating the 
regime below which this strong coupling description is 
appropriate. 

Our main finding, atleast for ^ 2/^1 < 1.5, as discussed 
below, is that the gapped spin liquid phase found in the 
slave rotor mean field theory in Section III may be poten¬ 
tially unstable to commensurate magnetic orders such as 
Neel or Tetrahedral order. On the other hand, the regime 
we identify as a gapless spin liquid with Fermi pockets 
in the slave rotor mean field theory appears to be domi¬ 
nantly unstable to various incommensurate spiral phases 
such as Spiral or Cone phases. Since the phase diagrams 
obtained in Fig. 2 and Fig. 3 involve different approxima¬ 
tions, neither one is fully reliable; further numerical work 
is necessary to definitely identify whether the spin liquid 
phases are indeed stable or give way to magnetic or other 
symmetry breaking orders. However, the true phase di¬ 
agram is expected to display phases shown in the phase 
diagrams in Fig. 2 and Fig. 3. 

We describe these various competing magnetic ground 
states below. In the Appendix, we include “Common 
origin” plotJ^ for these states for further insights. 

(i) Neel : The Neel state is collinear with spins pointing 

opposite to each other on the two sublattices of the brick- 
wall lattice, with (xa) = 0 on all triangular plaquettes. 
Its structure factor iS(q) = ^ ex¬ 

hibits a Bragg peak at q = (7r,7r). 

(ii) Tetrahedral : In this state, the spins point from the 
origin towards the four corners of a tetrahedron, and are 
tiled on the lattice as shown in Fig. This noncoplanar 
state has a uniform chirality (xa) = —1/6\/3 on each 
small-A. The tetrahedral state is a triple-q state, with 
the structure factor S{q) exhibiting Bragg peaks at q = 
±(7r/2, ±7r/2) and q = (0, tt). 

(iii) Cantellated Tetrahedral : The Cantellated Tetrahe¬ 
dral (CT) state descends from a parent tetrahedral state. 
It is reminiscent of the cuboctahedral state in frustrated 
kagome antiferromagnets.l^^M^ In the parent tetrahedral 
configuration, the spins marked A (or B^C^D) in Fig. 
form a honeycomb (or brickwall) lattice, with a larger 
unit cell. We can deform these ferromagnetically aligned 
spins by dividing the larger honeycomb lattice into two 
sublattices, and then allowing for a \/3 x \/3 canting on 
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FIG. 3: Top: Phase diagram of i/spin showing various classical 
magnetically ordered ground states at t//t = 10: Neel, Spi¬ 
ral, Tetrahedral, Cantellated Tetrahedral (CT), Cone-I, and 
Cone-IL Below the dashed (white) line, the slave rotor re¬ 
sult shows that we are in the Mott insulator, so the strong 
coupling expansion is expected to be valid. Stars indicate 
points where we numerically study the effect of thermal fluc¬ 
tuations. Bottom: Spin configurations in the Tetrahedral and 
Cone states (viewing the honeycomb lattice as a ‘brickwall’ 
for visual convenience). The indicated Triad states act as 
parent states for Cone-I and Cone-II, which are obtained by 
spiralling the Triads about the indicated hs axis. The CT 
descends from the Tetrahedral by a x canting pattern. 


angle 0 varies with ^ 2 , 0 ; at large ^ 2 /^ 1 , the spins on 
each large-A tend to form an orthonormal triad, with 
0 = tan“^(l/v/ 2 ). The Triad-I is a triple-q state; S{q) 
exhibits Bragg peaks at q = (0,±27r/3), q = ( 7 r,± 7 r/ 3 ), 
and q = (0,0). Let us now consider an axis hs which is 
perpendicular to the plane formed by h and any one of 
the spins. The Cone-I state is obtained by rotating all 
the spins about this axis by an angle ^(r) = q • r where 
q is an incommensurate spiral wavevector; moving along 
q, each spin of the triad again traces out a cone, and the 
Bragg peaks acquire a weak incommensuration. 

(v) Cone-II : The Cone-II state occurs for 0 > 7 r/ 3 , 
where the sin 30 term in Eq. changes sign. Here, 
the parent state is a Triad-II state in which the spins 
on sublattice-A are reordered, with = cos Oh + 
sin^ cos( 47 rj/ 3 )n^i + sin^ sin( 47 rj/ 3 )n^ 2 : while the cone 
on sublattice-B is simply flipped, so that = — cos Ohh- 
sin 0 cos((/:?+27rj/3)nxi -fsin 0 sin((/:?+27rj/3)n_L2, with the 
ground state energy being independent of ip. The Triad- 
II state has a nonzero staggered magnetization, with 
Bragg peaks in 5(q) at q = (0,±27r/3), q = ( 7 r,± 7 r/ 3 ), 
and q = ±( 7 r, 7 r). The Cone-II state is obtained by ro¬ 
tating all the spins of Triad-II about the hs = h axis by 
^(r) = q-r where q is an incommensurate spiral wavevec¬ 
tor, leading to weakly incommensurate Bragg peaks. 

(vi) Noncoplanar Spiral: At 0 = 0, the Hamiltonian 
L^spin is the frustrated Ji-J 2 Heisenberg antiferromagnet 
on the honeycomb lattice. The classical ground states 
of this model consists of coplanar incommensurate spi¬ 
rals for J 2 > Ji/ 6 .^ For 0 7 ^ 0, chiral terms in iLgpin 
cause spins on the A and B sublattices to cant in oppo¬ 
site directions away from the spiral plane, leading to a 
noncoplanar spiral with a uniform xa 7 ^ 0 on all big-As. 


V. THERMAL FLUCTUATIONS 


each sublattice; this splits each vertex of the tetrahedron 
into six vertices (forming a small hexagon) leading to the 
CT state with a 24-site unit cell. Here, the CT state ex¬ 
ists over an extremely narrow window between the Neel 
and Tetrahedral states. Third-neighbor • Sj terms, 
which favor the cuboctahedral state on the kagome lat¬ 
tice, could enhance the regime of CT order. 

(iv) Cone-I : This state descends from a Triad-I state 
in which spins on the large-A tend to form a triad on 
each sublattice, leading to a nonzero chirality. The or¬ 
ganization of the spins is shown in Fig. In Cone-I, 
which occurs for 0 < 7 r/ 3 , we start with triads on the 
two sublattices that are lined up with a common axis 
n, with spins making an angle 0 with h. Forming an or¬ 
thonormal basis {n, n_Li, ^± 2 }, we parameterize the three 
spins of the triad on sublattice-A, {j = 0 , 1 , 2 ), as 
= cos + sin^ cos( 27 rj/ 3 )nAi + sin^ sin( 27 rj/ 3 )n_L 2 , 
while the triad on sublattice-B is rotated by ip about n, 
with = cos0h^sm0cos{ip^27rj/3)h_\_i^sm0sm{iph- 
27rj/3)nA2- This state has a nonzero net magnetization. 
The ground state energy is independent of ip. The triad 


In 2D, the Mermin-Wagner theorenP^ precludes spin 
SU{2) symmetry breaking at T > 0. However, discrete 
orders associated with lattice symmetry breaking can sur¬ 
vive thermal fluctuations. 

The Neel and Tetrahedral orders break SU{2) symme¬ 
try, but their SU{2) invariant spin correlations, such as 
(S^ • Sj) and (S^ -87 x 8 /^), respect all lattice symmetries. 
These states at T > 0 can thus be smoothly connected to 
the high-T paramagnetic state. As shown in Fig. ga), for 
the Tetrahedral state at U/ti = 10, there is a crossover 
temperature ^ 0.05ti below which the spin chirality xa 
on small triangles becomes large, saturating to l/6^/3 as 
T ^ 0. The T > 0 regime here can thus be viewed as a 
elassieal CSL. By contrast, xa vanishes at T = 0 in the 
Neel state, but it exhibits a peak at nonzero temperature. 

Previous work has shown that at 0 = 0, the spiral 
ground states of the frustrated J 1 -J 2 Heisenberg anti¬ 
ferromagnet have nematic order associated with broken 
C 3 honeycomb lattice symm etry which survives up to 
a thermal Zs-clock transition.^ The noncoplanar spiral 
does not involve any additional symmetry breaking, so it 
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FIG. 4: (a) Temperature dependence of the spin chirality xa 
on small triangles for f//t = 10 and 0 = 7 r/ 3 , showing its low 
temperature saturation in the Tetrahedral phase — 1) 

and its vanishing in the T = 0 Neel phase (t 2 /ti = 0.4). We 
have plotted 5 x xa in the Neel state for clarity, (b) Log 
plot of finite size scaling of critical susceptibility for the Z 3 - 
clock ordering transition in the Cone-II state at Vjtx — 10, 
0 = 97 r/ 2 O, — 1.5 showing X 3 with 7 /?/= 1.80(8). 


undergoes a similar Zs-clock transition at 0 0 . 

In the Cone states, the choice of the common axis 
h represents a spontaneously broken SU{2) symmetry; 
long wavelength fluctuations will disorder h at arbitrar¬ 
ily small T > 0, restoring SU(2) symmetry. Using Monte 
Carlo simulations, we find that the nearest neighbor spin 
correlations (S^ • Sj) are modulated in the Cone states, 
leading to energy modulations on the bonds of the lattice. 
In the Cone-I state, these energy modulations resemble 
a \/3 X \/3 columnar valence-bond solid pattern, while in 
the Cone-II state the energy modulations on the bonds 
resemble a staggered valence-bond solid bond pattern. 
This leads to a Z 3 symmetry breaking in both states. 
We thus expect the Cone states to undergo a thermal 
Zs-clock transition into the high temperature paramag¬ 
netic state. In the Cone-II state at I7/ti = 10, (/) = 97 r/ 20 , 
^ 2 /U = 1-5, the computed peak susceptibility Xsi^) 
this transition is plotted in Fig.|^b); its finite size scaling 
shows Tc ~ O.lOti, with X 3 (U) with 7 /z^ = 1.80(8), 

consistent with the exact Zs-clock resulfP^ for the expo¬ 
nent ratio 7 /z^ = 26/15. 


VI. MOTT INSULATOR OF SPINOR BOSONS 
IN THE HALDANE HUBBARD MODEL 

Cold atom experiments can also study interacting 
bosons in the Haldane-Hubbard model. Previous work on 
spinless bosons has shown the emergence of chiral super¬ 
fluids and plaquette Mott insulators with loop currents.!^ 
For pseudospin- 1/2 bosons, the Haldane-Hubbard model 
is 


Here, we focus on the strong coupling limit 
Ulti,Ult 2 ^ 1. The spin Hamiltonian in the result¬ 
ing Mott insulator can be derived in a manner similar to 
the fermion case, although it is more tedious to calculate 
the terms to 0(t^/U^). We find 

s* • S,- + J2(<^) • Sj 

{ij) m) 

+ JJ2 N + ~ip N 

small—A big—A 

with Ji(0) = — 24^^ COS 0 and J2(0) = ~ 

6-^ cos 0—12^ cos30. Similar to the fermion case, the 
set of triangles includes small triangles that enclose flux 
= —05 big triangles within a hexagon which enclose 
a flux = 30, and big triangles defined around a site 
which enclose a flux ^a = —30. As expected, the two- 
spin exchanges are ferromagnetic at 0 = 0, and the chiral 
terms also change sign relative to the fermions. 

In order to localize bosons, it is well known that the 
required repulsive interactions are much stronger. Thus, 
t/U is much smaller in the Mott insulating phase of 
spinor bosons, so that the chiral terms are much weaker 
in magnitude compared to the fermion case. In addition, 
unlike the fermion case, both Ji and J 2 are ferromag¬ 
netic, which means these two-spin exchange couplings do 
not frustrate each other. Due to both these reasons, we 
find that ferromagnetic order almost completely domi¬ 
nates the magnetic phase diagram in the Mott insulator 
of spinor bosons. 

In the ferromagnetic Mott insulator, the spin of the 
bosons plays no role. Thus, our model should yield re¬ 
sults identical to that for spinless bosons. Indeed, assum¬ 
ing ferromagnetic order in our model, the chiral terms 
vanish. Nevertheless, we find that the bond energies dis¬ 
plays a flux dependence arising from the flux-dependence 
of the two-spin exchange couplings. This leads to a weak 
nonzero loop currents ^ /U‘^ around triangular plaque- 

ttes, consistent with weak loop currents predicted in the 
plaquette Mott insulator of spinless bosons using quite 
different approaches.!^ This result for spinor bosons is in 
striking contrast to the fermion case, where the only flux 
dependence arises from the chiral terms, so loop currents 
to 0{t^ IU‘^) only arise in fermion Mott insulators which 
have non-coplanar spins with a nonzero scalar chirality. 

VH. DISCUSSION 


= -iiE {blbj, + h.c.)-t2j2 {e^^^^*blb^. + h.c.) 

{^3)0- {{i 3 ))o- 

+ ( 3 ) 

iaa' 

For atoms such as ^^Rb, the background scattering 
lengths are nearly isotropic, so we set Uaa' = U. 


We have discussed spin- 1/2 fermions or bosons in the 
strongly correlated Haldane-Hubbard model. Our main 
finding is the rich variety of competing chiral magnetic 
orders with varying flux and hopping. Such orders, 
even if they are short-r anged , could be probed via mag¬ 
netic Bragg scatteringwhich has been shown to 
be a useful probe in atomic gases. For fermions, since 
dHspin/dcl) ^ Xa, measuring the energy change upon an 
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adiabatic change of flux is an indirect way to measure y a • 
This relation also ties the spin chirality to (weak) charge 
currents in the Mott insulator which could potentially 
be detected using quantum quenches More direct 
routes to detect the spin chirality of atoms are desirable. 
Doping these magnetically ordered states leads to non¬ 
trivial Chern bands and a nonzero charge Hall effect, sim¬ 
ilar to that in certain frustrated metallic magnet 
Quantum melting the Tetrahedral state, which is uni¬ 
form with a large chirality, may lead to a quantum CSL 
ground state with short range magnetic correlations at 
the same wavevectors. Further work is needed to clar¬ 
ify the energetic competition between CSLs and chiral 
magnetic orders. Finally, although the Mott insulating 


phase of spinor bosons is dominated by the ferromagnetic 
phase, lowering U might lead to superfluids with exotic 
magnetic orders; this will be explored elsewhere. 
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Appendix A: Common origin plots for the various magnetic orders 

For visualizing spin configurations, we employ common origin plotwhere all spins on the lattice are plotted with 
their tails at the center of a unit sphere and heads marked as dots on the surface of the sphere. We mark the spins 
on the two sublattices of the honeycomb (i.e., brickwall) lattice using different colors, red and black. Doing this for 
the Tetrahedral state, for instance, we see that all spins point along one of four directions, towards the vertices of 
a tetrahedron. For the Cantellated Tetrahedral, the vertices split into six points, three on each sublattice. For the 
Spiral state, we find that the spins on the two sublattices which support incommensurate spirals cant out of the plane 
of the coplanar spiral in opposite directions. The Cone-I state, which is obtained by spinning a triad state about a 
specific axis leads to one great circle and two small circles, while Cone-II leads to two parallel small circles. 





